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This paper is devoted to a derivation of the probability distribution of photon escape from a semi-infinite 
random medium, depending on the number of its interactions with macroscopic particles inside the medium. 
The consideration is limited to the case of highly developed multiple light scattering. The distribution function 
found facilitates the solution of both direct and inverse problems in light scattering media optics. 
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I. INTRODUCTION 


Statistical characteristics of a photon gas in random media 
is an important subject which has a number of applications 
[1-3]. We address the following question. We have a con- 
stant infinitely broad monochromatic photon flux incident in 
the direction specified by the unit vector mo (Ùo, Po) on the 
surface of a semi-infinite plane-parallel random isotropic me- 
dium. Here o is the zenith angle and @p is the azimuth 
angle of the incident light beam. It is assumed also that the 
energy of the photons is far from the absorption bands of 
substances contained in the random medium. We consider 
the angular distribution No(m, mgo) of photons emerging in 
the direction specified by the unit vector m(0,g) from an 
arbitrary point S on the surface of a layer. Here Ò is the 
zenith angle and ¢ is the azimuth angle of the emerging light 
beam. Clearly, due to the symmetry of the problem this dis- 
tribution does not depend on the choice of S on the surface of 
a random isotropic medium. Also we have for normal illu- 
mination mo: X=0, Xe L, where L is the plane containing the 
medium surface and x is an arbitrary unit vector in the plane 
L. 

The function No(m7,m9) can be obtained by solving Am- 
bartsumian’s nonlinear integral equation [4,5]. We are inter- 
ested, however, in the representation of No(m,my) as a sum 
of contributions p(m,7,,n) due to photons scattered n times 
in the medium [1-3]: 


co 


Noñ, iño) = D p(tii,tiig,n). (1) 


Photons scattered different numbers of times, of course, 
will have different path lengths in the scattering medium. 
The value of p(m,mgo,n), which is often called the photon 
weight, can be considered also in the framework of the Feyn- 
man path integral approach [6,7]. This approach, being very 
general, can be applied to any type of scattering medium. We 
will use here, however, the essential features of the medium 
under consideration, namely, its infinite extension in the 
space below the plane L and the absence of photon absorp- 
tion. These assumptions allow us to avoid path integral cal- 
culations. Clearly, the average number of scattering events in 
such an artificial medium is infinite as well. The primary 
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goal of this paper, therefore, is to derive the analytical rep- 
resentation for the weights p(m, 179 ,n). 


Il. THE PROBABILITY DISTRIBUTION FUNCTION 


Let us introduce the probabilities 
p(m,mo,n) 
f(n)= = (2) 
> pomon) 


n=1 


with the normalization condition 


X fn)=1. (3) 


n=1 


The value of f(n) can be interpreted as the probability for 
photons injected in the medium in the direction mo and scat- 
tered n times to emerge in the direction specified by the 
vector m. The condition (3) states that the total probability of 
photon escape is equal to 1. This is due to the assumed 
absence of absorption. 

To derive the function f(n) we will use the random walk 
theory [8]. This theory states that the probability of a particle 
appearing at a given place, time, and direction after a large 
number n of interactions is given by 


f(n)= Val mn expl- aln), (4) 


where the constant œ depends on the physical process under 
study. The only problem left is, therefore, to find the constant 
æ for our particular case. Clearly, it does not depend on the 
position of the point S on the surface of the medium. It also 
does not depend on time because we consider the steady 
case. So the only dependence left is due to local optical 
characteristics of the random medium and the vectors m, mo. 

To derive the parameter a, we consider now the case of an 
absorbing medium with the same scattering law in a single 
scattering event as for the nonabsorbing semi-infinite random 
medium in question. 

The probability of photon survival wọ in a single scatter- 


ing event differs from 1 for absorbing media. Then we have 
instead of Eq. (1) [9,10] 
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[ee] 


N(ñ,mo)= >, pñ, ñon). (5) 


n=1 


Clearly, it follows that N=No at wọ= 1. We can also use the 
distribution function f(n). Then we have for €=N/No 


oa fnot, (6) 


where for the sake of simplicity we omitted the indications 
of the dependencies of é and f on m, mo. We will consider 
now a random turbid medium with wọ™1. 

First of all we note that the expansion (6) is only slowly 
convergent for wọ™1. So we will use an expansion in the 
parameter £= 1— wo, which is the probability of photon ab- 
sorption, instead of the expansion in wọ in Eq. (6). Then we 
have from Eq. (6) 


Hà f(n)(1—B)" (7) 
or 


o0 X -] 
i= f| 1- gn OY) 


3 = = 
Bn(n ul 2) h m 
where we used the expansion 
a-p)"=> 1") (9) 
j=0 J 


with (n/j)=n!/j\(n—j)!. 
It follows from Eq. (8) that approximately 
Bn? Ben? 


*~exp(— Bn), (10) 


where 


nk= > f(n)n*, k=1,...,, 


exp(— Bn)= © f(n)exp(— Bn). (11) 


Note that we have assumed that n(n— 1)~n?, 
n(n—1)(n—2)~<n’,... in the derivation of Eq. (10). This is 
possible due to the large number of scattering events n. For 
the same reason we have 


TETE f Fonex- pnan, (12) 


where [see Eq. (3)] 
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[fwan (13) 


The integral (12) can be evaluated analytically. The answer is 


exp(— Bn)=exp(—2V\aß), (14) 


where we used Eq. (4). Thus we have, taking into account 
Eq. (10), 


N(m, Mmo) =No(m,mo)exp(—2 VaB) (15) 
or, as B—0, 
N (iii, tig) =No( ti, mig) (1-2 Va). (16) 


Now we have an opportunity to find the value of a, com- 
paring Eq. (16) with the exact solution of the radiative trans- 
fer equation at small £ [3]: 


N(m,mo)=No(m,mo)[1—yu(m, mo) ], (17) 


where y=448/3(1 — g) and g=3JSjp(@)sin Ocos 6d6 is the 
asymmetry parameter [3], p( 0) is the probability of a photon 
scattering in the direction specified by the scattering angle 6, 
and 


u(m,mo)=[Ko(“)Ko( Mo) V/[Ro(u.M0,¥)], (18) 


where [1] 


plz 
Kow=5 | Ra( u, Ho) (M+ Mo)d po (19) 


is the escape function and 


1 2m 
Reno) === | RS (u, Mo. Wd (20) 


is the azimuthally averaged reflection function R?( H, Ho Y) 
[3]. Here u=cos 3, uo=cos %, P=E- go. Clearly we 
have y= V1- (ño 8, u=V1-(m- xX). 

It follows from Eqs. (16) and (17) that 


a=(4u7)/[3(1—g)], (21) 


which is the result we tried to establish from the very begin- 
ning. Finally, we have from Eqs. (4) and (21) 


2u exp[ —4u7/3n(1—g)] 
n? 3 m(1—g) 


where u is given by Eq. (21). We note here the importance of 
the viewing function u, which combines all angular depen- 
dencies. 

Our Eq. (22) transforms to the similar equations derived 
in [9,10] if one uses the expansion of the exponent in a 
power series in the value of n~!. Only the first term of such 
an expansion was explicitly derived in [9,10]. 

It follows from Eq. (22) that f’(n)=0 at n=n,yyx=2a/3 
and 


f(n)= > (22) 
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FIG. 1. The dependence f(n) for various values of g at u=1. 


1 3 3/2 3 
J max=S (nmax) = alm | ;| ox ~~ ;| > (23) 


where a is given by Eq. (21). We also have [see Eq. (21)] 
Nmax= (8u")/[9(1 —g)]. (24) 


We see that nmax increases with g and fmax is linearly pro- 
portional to |—g. 
Finally, it follows from Eqs. (2) and (22) 


P(H, Hon) 


_ 2Ko(M)Ko( Ho)expl — 4u? (u, mo, W)/3n(1— g)] 
n>? /3a(1—-g) í 


(25) 


where we used the normalization condition No(m,mo) 
=Ro(u, uo, Y). The function Kọ(u) can be approximated 
by [11,12] the expression Ko()=3(1+2,) for arbitrary 
random media with discrete particles. The accuracy of this 
formula is better than 2% at w>0.2 [11,12]. Simple approxi- 
mations for Ro(, uo, Y) are derived in [1,5,11,12]. 

We present results of calculations with Eq. (22) at u= 1 in 
Fig. 1. It follows that the maximum of the probability f 
increases linearly with increasing 1 —g, which is in corre- 
spondence with Eq. (23). Clearly, we obtain that f=0 at g 
=|. This means that photons do not have a chance to escape 
from the medium. They only propagate along straight lines 
(no scattering). 

The maximum shifts to larger numbers of scatterings n for 
larger values of g, which is in correspondence with Eq. (24). 
The influence of the parameter s=u 7 on the curves f(n) is 
similar to that of the parameter 1 — g. The value of nmax only 
slowly changes with g at g<0.7. However, it increases rap- 
idly at g>0.9. 
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HI. APPLICATIONS 
The main result, given by Eq. (22), can be generalized to 
other situations and measurement setups. For instance, if one 


is interested in functions f(n) for total photon numbers, in- 
tegrated over the escape angle, it is possible to obtain instead 
of Eq. (22), following the same line of reasoning, 


Kn- 2K o(Mo)expl — 4K o(Mo)/3n(1 —g)] 
n*?3 «(1 —g) 


Using similar arguments we obtain for the case of inte- 
gration on both the incident and escaped photon directions: 


(26) 


= 2 exp[ —4/3n(1—g)] 
f=, == 
n°3 (1-8) 
This corresponds to the case of diffuse illumination and dif- 


fuse reflectance measurements. Clearly, we have for the dif- 
fuse reflection coefficient r( 09) 


(27) 


r( Bo) > f(n)@o, (28) 
which gives us after transfer to the continuous basis 
r(9o)= [ Totan. (29) 
This formula can also be presented in the following form: 
r(0)= | “Flmdexp(pnan, G0) 


where p=In(1/@9)~1—@ =P as wo— 1. It follows from 
Eqs. (30) and (26) that 
r( Bo) =expl —yKo(Ho)], (31) 


where y =4 y[In(1/a 9) ]/3(1 —g). Also, we have for the total 
reflectivity 


> fas (32) 


or, following the same steps as in the derivation of Eq. (31), 
r=exp(—y). All cases considered here correspond to weak 
absorption and, therefore, we can use the limiting value of y 
for wy close to 1: y =4/(- @,)/3(1—g). The expres- 
sions for r(%) and r derived here have been known for a 
long time [11]. They were obtained, however, using a differ- 
ent approach. Our derivations allow us to make clearer their 
physical basis. 
Also, using Eqs. (22) and (6) we conclude that 


_ Ro 
¿= Ro =exp(—uy). (33) 


This formula is also the well known result of radiative trans- 
fer theory [11,12]. We see that the derived function (22) ap- 
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pears to be a key point for the derivation of many important 
relations in scattering media optics in a simple and straight- 
forward manner. 

It can also be used to establish the temporal statistics of 
photons emerging from a random light scattering layer. For 
instance, accounting for the fact that n=vt/L, where t is 
time needed for a photon travel the distance nL with the 
group speed v, we may obtain from Eq. (22) the distribution 
of photons according to arrival times f or paths nL. Here L 
=o _!, where ø is the extinction coefficient of the random 
medium. The value of L is called the photon free path length. 
Taking this into account, we obtain from Eq. (22) 


E t 3/2 


Galan 


Clearly the functions f(n) and Fn) transform to 


-o 
3(1—g)vt 
(35) 


| 4Lu? 
P 3(1—g)ut 


U 


fle)=2u| | (34) 


t 3/2 


ee 


L 3/2 
Fl)=2Kol40)( E 


and 


B L\32 p2 
a| | (36) 


4L 
sl r 


For instance, we have from Eq. (35) at normal incidence 


_ A a BL 
f(t)= y[B/7(1 ail <| ox oan} (37) 


where 


U 


B= K}(1). (38) 


The value of B is approximately equal to 2.2. Equation 
(37) has been derived earlier [7,13]. However, the value of B 
in [13] is equal to 0.75. We believe that B, given in Eq. (38), 
is closer to the exact result. Note that the value of B is equal 
to 2.19 in [7], which is close to the result given by Eq. (38). 

The distribution on path lengths is obtained by substitu- 
tion of t in Eqs. (34)-(36) by s/v, where s=nL is the total 
distance traveled by a photon after n scatterings. 
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The probability distribution function f(n) can be used 
also to find the statistical moments 


nk= [ faan. (39) 


0 


It can be done analytically; namely, we have 
k-1 
F =i ð exp(—2vVap) 
n*=(—1) a4 — 
ap Vp 


where a is given by Eq. (21) and p=In(1/w9). In particular, 
we have 


| » (40) 


n= alp exp(—2Vap). (41) 


Thus, one can obtain for the average number of scatterings 
involved in forming the observed absorption line [14] 


in| [nemo] / | [ferret (42) 


the following simple relation: 
(N)=(2u)/¥, (43) 


where y= V3 (1—g)(1—@ 9) is the diffusion exponent of ra- 
diative transfer theory [3]. The average distance traveled by a 
photon having (N) scatterings is given by /=(N)L. This dis- 
tance varies with the observation geometry due to the pres- 
ence of the viewing function u in Eq. (43). Thus, we con- 
clude that the strength of the absorption line in the scattering 
atmosphere will also depend on the viewing geometry. 


IV. CONCLUSION 


In conclusion, we derived here the probability distribution 
function (22). It describes the photon migration from the 
direction mp to the direction m after n interactions with scat- 
terers. This function might be of importance for a wide range 
of applications. Some of them are outlined above. Interest- 
ingly enough, the geometry of the observation enters Eq. 
(22) as a single number u. This gives significance to the 
function u, given by Eq. (18), for radiative transfer theory. 
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